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We present an efficient quantum algorithm for preparing a pure state on a quantum computer, 
where the quantum state corresponds to that of a molecular system with a given number m of 
electrons occupying a given number n of spin orbitals. Each spin orbital is mapped to a qubit: the 
states |1} and |0} of the qubit represent, respectively, whether the spin orbital is occupied by an 
electron or not. To prepare a general state in the full Hilbert space of n qubits, which is of dimension 
2", 0(2") controUed-NOT gates are needed, i.e., the number of gates scales exponentially with the 
number of qubits. We make use of the fact that the state to be prepared lies in a smaller Hilbert 
space, and we find an algorithm that requires at most 0(2'"^^n'"/m!) gates, i.e., scales polynomially 
with the number of qubits n, provided n S> m. The algorithm is simulated numerically for the cases 
of the hydrogen molecule and the water molecule. The numerical simulations show that when 
additional symmetries of the system are considered, the number of gates to prepare the state can be 
drastically reduced, in the examples considered in this paper, by several orders of magnitude, from 
the above estimate. 
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I. INTRODUCTION 

Simulating quantum systems on a classical computer 
is a hard problem. The size of the Hilbert space of the 
simulated system increases exponentially with the sys- 
tem size. For example, in quantum chemistry, the full 
configuration interaction (FCI) method diagonalizes the 
molecular Hamiltonian to provide solutions to the elec- 
tronic structure problem. The resource requirement for 
performing the FCI calculation scales exponentially with 
the size of the system Therefore, it is restricted to the 
treatment of small diatomic and triatomic systems 0- 
Feynman 0] observed that simulating a quantum system 
might be more efficient on a quantum computer than on 
a classical computer. Further work has born out of this 
early suggestion 0, i, i, 0, i, i M EH, [H, [11, El • 

The implementation of a quantum simulation algo- 
rithm requires a mapping from the system wave func- 
tion to the state of the qubits. One possible mapping is 
through the Jordan- Wigner transformation (JWT) [15j . 
Each spin orbital is mapped to a qubit: the states |1) 
and |0) of the qubit represent, respectively, whether the 
spin orbital is occupied by an electron or not. The simu- 
lated system, and therefore the quantum computer, can 
be in a quantum superposition of different configura- 
tions. In quantum chemistry, this quantum superposition 
is known as the configuration state function (CSF). There 
are other possible mappings between the simulated sys- 
tem and the quantum computer [l,E3l- -f'o'" example, one 
can make use of the fact that the Hilbert space of inter- 
est, i.e., one with m electrons occupying n spin orbitals, 



is of dimension 



and use the smallest number n' 



of qubits that satisfies the condition 2" > 



The 



large enough to describe any quantum state of the sim- 
ulated system. However, since these compact mappings 
complicate the implementation of the simulation, we only 
consider the simple mapping explained above. 

Preparing a general pure state in the full Hilbert space 
of 71 qubits, which is of dimension 2", has been studied by 
several groups. Shende and Markov [l^ and Mottonen 
et al. [17| showed that preparing a generic n-qubit pure 
state from |0)®" requires 0(2") controUed-NOT (CNOT) 
gates, which gives an exponential scaling for the num- 
ber of CNOT gates; Bergholm et. al. gave an up- 
per bound for the number of gates required for trans- 
forming an arbitrary state \a) to an arbitrary state \b): 
(2"+i-2n-2) CNOT gates and (2"+i-?i-2) one-qubit 
gates. The gate count is halved if \a) or |6) is one of the 
basis configuration states. Another quantum algorithm 
for the preparation of an arbitrary pure state with fi- 
delity arbitrarily close to one was suggested by Soklakov 
and Schack [19[, which is based on Grover's quantum 
search algorithm. 

In this paper, we study the preparation of a pure state 
in the second-quantized representation that represents 
the quantum state of m electrons distributed among n 
spin orbitals. The state to be prepared lies in the combi- 

n 



natorial space of dimension 



which is a subspace 



Hilbert space of the quantum computer would then be 



of the full Hilbert space of n qubits. Ortiz et al. [20[ 
studied a similar problem, and they gave an algorithm 
for preparing a state that is composed of N configura- 
tions, where TV is a finite and small number. Their algo- 
rithm scales as 0{N'^v?), where n is the number of the 
qubits. In this paper we present an efficient recursive 
algorithm that gives specific quantum circuit for prepar- 
ing a pure state, with polynomial scaling of the number 
of CNOT gates in terms of the number of the qubits. 
We have numerically simulated the state preparation al- 



gorithm to prepare the electronic states of the hydrogen 
and the water molecules. The results show that the num- 
ber of CNOT gates can be reduced by up to 3 orders of 
magnitude from the upper bound we derive. 

The structure of this work is as follows. In Sec. |TI] 
we discuss mapping the Fock space of the system onto 
the Hilbert space of the qubits. In Sec. IIIII we present 
a recursive algorithm for state preparation. Preparing a 
general pure state in the one-electron systems is discussed 
in detail. Then based on the results for the one-electron 
system, a recursive approach for preparing a general state 
of an TO-electron system is presented. In Sec. IV, we dis- 
cuss the simplification of the two-fold controlled unitary 
operations in the algorithm. In Sec. V, we analyze the 
scaling of the algorithm. In Sec. VI, we apply the algo- 
rithm to prepare the electronic states of the hydrogen and 
the water molecules. We close with a conclusion section. 



II. FOCK SPACE OF MOLECULAR SYSTEMS 
AND THE JORDAN- WIGNER 
TRANSFORMATION 

The wave function of a m-electron, n-spin-orbital sys- 
tem is a linear combination of Slater determinants. In 
the formalism of second quantization it is written as: 

l*>= E K,■■■,^Mr■■■{aiy-\n), (1) 

ii, ■■■ , irze{o,i} 

il + ■■■ -f- in—m 

where jfi) represents the vacuum state with no elec- 
trons. Each set {zi, • • • , i„} represents an electron con- 
figuration of the system, and these states are used as 
the basis vectors of the configuration space of the sys- 
tem. The fermion creation and annihilation operators 
and Oj satisfy the canonical anticommutation relations 
{ai,aj} = 0;{ai,aj} = (5y. (a^) creates (annihilates) 
a fermion on the jth spin orbital. 

The fermion algebra is isomorphic to the standard 
quantum computing (QC) model (or Pauli) algebra. The 
isomorphism is established through the Jordan- Wigner 
transformation [l^ . The JWT maps a fermion state to a 
one-dimensional standard QC state and vice versa. The 
creation and annihilation operators for the fermion state 
are mapped to the Pauli operators through the JWT, 

«.-[n£^(-cT^)]ai, (2) 

al-[n£^(-^z)]^;, (3) 

where cr^ is the Pauli matrix defined as: 

The operators and tT_|_ are the Pauli lowering and 
raising operators defined as: 
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= ^(O'x + i(Ty) = ) " ^^"^ 

They satisfy the following conditions: 

a+|0)=0, a_|0) = |l), (7) 

and 

a+|l) = |0), (T-|1)=0, (8) 

where 

|0)=(J), Il)=(?)- (9) 

III. RECURSIVE ALGORITHM FOR STATE 
PREPARATION 

State preparation means starting from a standard ini- 
tial state on a quantum computer and, by applying a 
unitary operation, transforming it to a given target state. 
In this paper, the target state in general is an entangled 
state mapped from a state in the Fock space of the m- 
electron, n-spin-orbital system. The unitary operation 
can be decomposed into a sequence of elementary quan- 
tum gates. 

Our approach for state preparation is as follows: in- 
stead of starting from the initial state and transforming 
it to the target state, we start from the target state and 
transform it to the initial state. We can obtain the gate 
sequence for preparing the target state by inverting the 
quantum circuit sequence, since any unitary operation is 
reversible. A recursive procedure is applied in our "re- 
verse engineering" approach. We first solve the problem 
of state preparation for one-electron system, then based 
on this, the general state preparation for m-electron sys- 
tem is solved. 

We set the initial state of an n-qubit system to be 10)*^" 
and define the target state as: 

N 

|*t) = E^»I^»)' (10) 

i=l 

where N is the dimension of the configuration space and 

N 
i=l 

where ji/'i) is a single configuration that describes a dis- 
tribution of m electrons on n spin orbitals. The states 
{iV'i)} are the basis vectors of the configuration space. 
For example, the state ^|001) + ^|010) + ^|100) is 
composed of three configurations in the state space of a 
one-electron, three-spin-orbital system. 
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FIG. 1: Quantum circuit for the controlled- iif gate. 



A. Gate library 



FIG. 2; Quantum circuit for transforming a Bell-type state 
1/3) = a|10) + 6|01)) to the initial state |00). So, |00> = ■ 
CNOT-X^ ■ 1/3). 



The gate library used in our approach contains the 
CNOT gate and the single qubit gates C. Here C repre- 
sents all unitary single qubit operators in SU(2). It can 
be written in the form: 



C = 



i;|l)(0|-t-*|0)(l|+«*|l)(l| = 



V u 



where u and v are complex numbers and 



\u? + \v? = \. 



(12) 



(13) 



The Hadamard gate H and the NOT gate X are included 
in C. We define another class of operators iJ, which is a 
sub-class of the single-qubit operator C, such that: 



H = C'^XC = 



U V + UV U — V 



-{u*V + UV*) 



(14) 



for some C. We call H the generalized Hadamard gate. 
It has the same property as the Hadamard gate, HH = I, 
and acts in a similar way as the Hadamard gate. Just 
like the Hadamard gate transforms a superposition state 
^ (|0) + 1 1) ) of a single qubit to state 1 0) and ^ ( 1 0) - 1 1) ) 

to state |1), the H gate transforms a given superposition 
state a|0) -I- &|1) of a single qubit to state |0) or |1). 

A considerable effort has been made to synthesize 
two-qubit circuits using CNOT gates and single-nubit 
gates_El Emm and quantum logic circuits [iH, S 
\2TI [28^ . It has been shown that to implement a typi- 
cal two-qubit operator, three CNOT gates are needed. 
In this paper we only use two kinds of two-qubit gates: 
the CNOT gate and the controlled-i7 gate, C-H, which 
as shown in Fig. 1, contains only one CNOT gate. In 
our algorithm, the H and C-H gates are used to address 
the configuration coefficients and reduce the number of 
the configurations that span the target state |5't) until 
it reaches the initial state jO)**". An example for trans- 
forming a three qubit state to the initial state is shown 
below. 

To transform the target state to the initial state, we 
first disentangle the target state, then reduce the num- 
ber of configurations in the target state by applying H 
and C-H gates. A simple example for a Bell-type state 
1/3) = a|10) -I- 6|01) is shown as follows: to transform 
this state to the initial state |00), first we apply the 
X^ gate to the first qubit (the number "1" on the su- 
perscript indicates that the operation is applied to the 
first qubit). We thus obtain X^\l3) = a\00) + b\ll). 



We then disentangle this state by applying a CNOT 
gate: CNOT-^'^ (where on the superscript, the first num- 
ber represents the control qubit and the second num- 
ber represents the target qubit). Therefore we now ob- 
tain CNOTi'2 . x^\(3) = a|00) + b\10) = (a|0) + 5|1))|0). 
Finally by applying a generalized Hadamard gate H^ 
to the first qubit, we obtain the initial state |00). So, 
|00) = H^ ■ CN0Ti^2 . x^\P). Thus the inverse gate se- 
quence will transform the initial state |00) to the target 
state 1/3). The quantum circuit for this process is shown 
in Fig. 2. 

The following is an example that shows how the con- 
figuration coefficients can be addressed and the number 
of the configurations in the target state can be reduced 
through the application of the generalized Hadamard 
gate H and C-H gate. To transform the state |5't) ~ 
-i=|001) + ^|010) + ^|100), to the state |000) = \0)®^, 
the gate sequence is as follows: 



V2' 



■V3 

CNOT^-\ |l)(i^|ll) + i^|00)) + i= 



V3 Ve y/2 

CNOT'^\ |l)(_L|io) + -L|oo)) + -L|o)« 

C-H''\, -^|100) + -i=|0)«3 



V2' 
H^ |0)^3^ 



(15) 



where in the C-H operation, the H gate is determined 
by solving the following equation: 




2uv u 



-2uv 



(16) 

We set u and v in the unitary matrix C to be real numbers 
since the coefficients in the target state are real. Then 
one can obtain the unitary matrix C by solving the above 

equation; one obtains u ~ i(3\/2 + 2\/3)v3 — \/6 and 



V ~ —yl/2 — l/\/6. The quantum circuit for this pro- 
cedure is shown in Fig. 3. 
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FIG. 3: Quantum circuit for transforming the three-qubit 
state |*T> = ^|001> + ^|100) + ^|010) to the state |0)®^ = 
1000). 



B. One-electron system 

Let 1)) denote the state wave function of the 

one-electron, n-spin-orbital system. If we consider aU 
possible distributions of the electron on the n spin or- 
bitals, the dimension of the configuration space will be 
n 



1 



— n. We call this the complete configuration space 



of the one-electron system. The target state can be fac- 
torized as follows: 

|*T(r^,l)) -co|0)|vl/T(n-l,l))+ci|l)|0)®("-i), (17) 

where |4't("-— 1, 1)) is a state wave function for one elec- 
tron distributed on (n — 1) spin orbitals. The configura- 
tion coefficients in state — 1,1)) are normalized 
to 1. Also, Co and Ci satisfy |cop + |cip = 1. The 
relative phase between two states can be addressed by 
some unitary operations. We define an unitary operator 
Q{n — 1, 1), such that 



Q(n-l,l)|vl/T(n-l,l)) = |0)«5("-i). 



(18) 



The procedure for transforming the target state to the 
state jO)*^" can be formulated as follows: 

\^T{n, 1)) coll) \^T{n - 1, 1)) + ci|0)®" 

C-g(n- l,l){i'(2,-.n)} co|l)|0)®("~^) +ci|0)®" 

> 

|0)®", (19) 



where C-Q{n — 1, l){i'(2,- - ,«)} jg ^ controlled unitary op- 
eration with the first qubit as the control qubit, and 
qubits 2, • • • , n are the target qubit. Here, the H gate op- 

Cl 

. Co 

Q ^ = |0). Thus the gate H used here is determined 
by solving the following equation: 



erates on the vector 



and transforms it to the state 



H 



Z) (o 

(20) 

The quantum circuit for this procedure is shown in Fig. 4. 



U*V + UV* U* — V* 

— f ^ —{u*v-^uv*) 
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FIG. 4: Quantum circuit for transforming a general state of 
the one-electron system |^T(?^, 1)) to the initial state |0)®". 



The unitary operator Q{n — 1,1) that transforms the 
state |*t("- - 1,1)) to the state |0)«'("-^' can be fac- 
torized similarly through the factorization of the state 

|*T(n-l,l)), 

|vI/T(n-l,l)) =c;,|0)|vI'T(n-2,l))+c;|l)|0)^("-2)^ (21) 



where |co| 



|2 



1 and \'i>T{n — 2, 1)) is normalized to 



1. An H' gate acts on the vector 



and transforms it 



to the state 



|0). The -ff gate can be determined 



by solving an equation that is similar to Eq. (20). We 
now define a unitary operator Q{n — 2, 1), such that 

Qin - 2, l)|*T(n - 2, 1)) = |0)«("-2)^ (22) 

The decomposition procedure above is repeated until 
the target state reaches a Bell-type state |\I't(2,1)) = 
Co|01) -I- c"|10), whose disentanglement procedure has 
been explained above (Fig. 2). The quantum circuit 
for this procedure of decomposing the unitary operator 
Q{n — 1, 1) is shown in Fig. 5. 



e(«-u)= 



- X 


1 
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e(n-2,l) 









H 



FIG. 5: Quantum circuit for the decomposition of the unitary 
operator Q{n ~ 1, 1). 

From the derivation above, we can see that the prepa- 
ration of the state |5't('t-, 1)) is reduced to the prepa- 
ration of the state |^t(2,1)), which is solved above as 
shown in Fig. 2. The Toffoli gate, which is also known as 
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the "controUcd-controlled-not" gate, appears when im- 
plementing this recursive procedure. We will discuss the 
simplification of the TofFoli gates into CNOT gates in Sec. 
IV. The scaling of the algorithm in terms of the number 
of CNOT and single-qubit gates is analyzed in Sec. V. 

States in the incomplete configuration space of the one- 
electron system are easier to prepare than those in the 
complete configuration space, since the former can be 
reduced to states in a complete configuration space of 
lower dimension. 



C. m-electron system 

For the m-electron system, it becomes more difficult 
to prepare the target state than for the one-electron sys- 
tem. We apply a recursive approach to solve this problem 
based on the following facts: 

(1) Any state of the one-electron system can be pre- 
pared [as shown in Sec. HI. B.] 

(2) An arbitrary state \'^t{p,P — 1)) can be prepared. 
[By applying X'^p gates, the state is transformed to 
\'^t{p, 1)), which is a one-electron state.] 

We now define a unitary operator U{n, m) such that: 



U{n,m)\'^Tin,m)) = |0) 



(gin 



(23) 



A quantum circuit for U {n, m) can be obtained through 
the decomposition of the state \'9T{n,m)) as follows. 

Any arbitrary target state [Eq. (10)] can be rewritten 
in the following way: 



The relative phase between two configurations can be in- 
cluded in these unitary operators. Then the target state 
can be transformed to the initial state jO)**" in the fol- 
lowing way: 

|4'T(n,m)) X\ co|l)|^r(n-l,m)) 
+Ci|0)|*T(«-l,m-l)) 

U{n- l,m - 1) 

> 

co|l)C/(n- l,m- 1)|*t(?i- 1,to)> +ci|0)®" 
C-C/(?i-l,m-l)-i^'-^''-'"» 



Co|l)|*T(?i-l,m))+ci|0)«" 
C-C/(n- l,m){i'(2'-'«)} 



co|l)|0)®("-i) +ci|0)®" 

, |0)®". (29) 

The quantum circuit for this procedure is shown in Fig. 6. 



U{n,m) = 
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-U[n-\,m-i)- 


U~\n-l,m-])- 


- u(n-l,m, 













III 



FIG. 6: Quantum circuit for transforming a general state 
\'^T{n, m)) of the m-electron system to the initial state jO)®". 



No Ni 

i=i j=i 

where the state is divided in two parts: one of them lies in 
the subspace where the first qubit is in state |0) (we call 
this the jo) subspace); the other part lies in the subspace 
where the first qubit is in state |1) (we call this the |1) 
subspace). Also, A^o is the dimension of the |0) subspace 
and A^i is the dimension of the |1) subspace. The target 
state can be factorized as follows: 

\-^T{n,m)) = co|0)|*T(n-l,m))+ci|l)|*T(7i-l,TO-l)). 

(25) 

The coefficients cq and ci are calculated as: 

|co|' = El^o,.|^ |ci|' = EI^-mI'' (26) 

1=1 3=1 

where jcoj^ -I- jcij^ = 1. We now define unitary operators 
U(n — 1, to) and U{n — 1, to — 1) such that 

U{n ~ l,m)\^T{n - 1, to)) = |0)®("-i), (27) 
C/(n- 1,TO- 1)|*T(?^- 1)) = |0)®("-i'. (28) 



The unitary operators C/(n — 1, to) and C/(n — 1, m — 1) 
that transform the states Y^riji — 1,to)) and \'^T{n — 
l,m — 1)) to the state jO)'^^""^-' can be decomposed 
through the decomposition of the states \^T{n — l,m)) 
and |^'t('^ ~ l,m— 1)). For example, for the operator 
U{n — 1, TO — 1), we have: 

\-^T{n-l,m-l)) = <^o\0)\'^T{n-2,m-l)) 

+c'i|l)|*T(ri-2,m-2))(30) 

Define unitary operators U{n — 2, to— 1) and U{n — 2, to — 
2) such that 

C/(n-2,TO- l)|*T(n-2,TO- 1)) = |0)®("-2)^ (31) 

C/(n-2,TO-2)|*T(n-2,TO-2)) = |0)^("-2). (32) 

The quantum circuit for U{n—l,m — l) is similar to the 
circuit in Fig. 6, where n and m in the unitary operators 
on the bottom line of the circuit arc replaced by (n — 1) 
and (to — 1), respectively. 

One can repeat this procedure until reaches a one- 
electron state |^'t(p, 1)), where 2 < p < n, or a 
I^t(P:P— 1)) state, which is solved as we discussed above. 
In each iteration, the unitary operators can be decom- 
posed into the circuits that have the same structure as 
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shown in Fig. 6. Such that the problem of transforming 
a general m-electron, n-spin-orbital state |5't(?t,, m)) to 
the initial state 10)*^" can be solved. This algorithm also 
works for systems of different numbers of the electrons. 
The simplification of the twofold controlled unitary op- 
erations and the scaling of the algorithm are discussed in 
Sees. IV and V. 



IV. SIMPLIFICATION OF THE TWO-FOLD 
CONTROLLED GATES 

Because controlled gates appear inside the controlled 
operations in Sec. IE, the twofold controlled gates ap- 
pear in the implementation of the recursive procedure 
for state preparation. For the simplest twofold controlled 
gates, the Toffoli gate [1^, which is also known as the 
"controUed-controUed-not" gate, in general, six CNOT 
gates are needed for implementing a three-qubit Toffoli 
gate [soj . However, if the input state for the quantum 
circuit is known, then the operation of a given unitary 
operator on all orthogonal states is immaterial. In such 
cases the unitary operator is said to be incompletely spec- 
ified p^ . Thus, there exists more than one quantum 
circuit for performing such an operation. So, one can se- 
lect a simpler circuit in order to reduce the number of 
gates. Based on this, we will now simplify the two-fold 
controlled gates in the algorithm. 

In our algorithm, the problem of preparing an n-qubit 



state is reduced to the problem of preparing some sim- 
pler states [i.e., for one-electron system, the problem is 
reduced to the preparation of a Bell-type state; for a 
general m-clectron system, the problem is reduced to the 
preparation of states \'^t{p, 1)) and |^'t(p,P — I))]- We 
will show that the twofold controlled gates that appear 
in each step of this recursive procedure can be simplified 
to onefold controlled gates by just turning off the opera- 
tion of the first control qubit. The recursive approach for 
constructing the quantum circuit leads to another advan- 
tage: for the controlled operations, it provides a way that 
makes the shortest distance between the control qubit 
and the target qubit since the controlled operations are 
simplified in each step while expanded to one more qubit. 

For the one-electron system, plugging the circuit for 
Q{n — 1, 1) shown in Fig. 5 into Fig. 4, the quantum cir- 
cuit for transforming the target state |\I't(?t-, 1)) to the 
state jo)**" is shown in the top part of Fig. 7. The quan- 
tum circuit has the following structure: from right to left, 
the input state is 10)**"; an H gate operates on the first 
qubit and followed by a C-H gate. After the application 
of these two gates, the basis vectors that span the state 
on the two control qubits are {|00), |10), The basis 

vector 1 01) does not appear, such that the two- fold con- 
trolled gates following these operations are incompletely 
specified. From right to left, taking jO)*^" as the input 
state, the basis vectors that appear in the state of the 
qubits change as follows: 



~ r io\io\«("-i) 1 .,12 lo>lo>|o>«("-^M ,^^,,3 , |o)|o)|o)«("-2) 



|l)|l)|0)®("-2) J y |l)|l)Q(n-2,l)-i|0)®("-2) 

J331 



The twofold controlled unitary operator C- 
(5~-'^(n — 2, l)'^'^"*''^' is incompletely specified. 

Here the use of C-Q-i(n - 2, l)'^'^^'^^' " instead of 
C-Q{n - 2, " is because we let the gates 
operate to the right. As a result, it can be simplified 
to C-(5~^(?i — 2, 1)^'^^' This simplification can be 

repeated in each recursive step until the state reaches 
a Bell-type state |4't(2,1)). The simplified quantum 
circuit is shown in the bottom part of Fig. 7. 

For the m-clcctron system, each of the three unitary 
operators in the bottom line of the quantum circuit in 
Fig. 6 can be decomposed into the circuit that has the 
same structure as that of U{n,m). The circuit in Fig. 6 
is divided in three parts. The two-fold controlled unitary 
operator in the controlled operations in parts 11 and /// 
will appear when they are expanded to one more qubit. 
Taking the state 10)*^" as the input state, from right 
to left; following the operation of the circuit, the basis 



vectors that appear in the state of the qubits change in 
the same way as that of the one-electron system. These 
twofold controlled unitary operators in parts E and IH of 
Fig. 6 are incompletely specified. For quantum circuit 
U{n — 1, m — 1) that acts on qubits 2, • • • , n, in part I of 
Fig. 6, the twofold controlled unitary operators may ap- 
pear when the circuit is decomposed one more step. From 
left to right in Fig. 6, after the circuit in part I operates 
on the target state |^'t('^,w)), the target state is trans- 
formed to ci|l)C/(n- l,m- l)|«'T(n- l,m)) -hco|0)®". 
This intermediate state can be divided in two branches: 
the |0) branch, lO)"^", and the |1) branch, \l)U{n-l,m- 
l)\'i'T{n — l,m)). For the |1) branch, the operation of 
U{n — 1, m — 1) that acts on qubits 2, • • • , n, will be can- 
celled with the controlled operation C-U~^{n — 1, to — 1) 
in part I. For the |0) branch of the state, the state 
on qubits 2, • • • , n, is lO)^*^""^-*, and it acts as the in- 
put state for the circuit U{n — 1,to — 1) from right to 
left. Then the unitary operator U{n — 1, m — 1) has the 
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H 



G(«-2,i) 



v®{»-2) 



H 



Q{n-2,l) 



in-2] 



FIG. 7: Top quantum circuit for the recursive procedure to 
transform the state \^T{n, 1)) to the initial state |0)®". It has 
two one-fold control operations and one two-fold operation in 
the middle. The bottom circuit shows its simplified version, 
without two-fold operations. 



same structure and input state as that of U{n, m) we dis- 
cussed above. The twofold controlled unitary operators 
inC/(n— 1,TO— 1) (parti of Fig. 6) are also incompletely 
specified. 

According to the analysis above, the unitary operator 
U{n — 1, m — 1) has the same structure as the circuit for 
U{n,m), with state jO)®*""^' as the input state. The 
twofold controlled unitary operators in the controlled 
unitary operators in [/(n — 1, to — 1) and U{n, m) are in- 
completely specified. A general form of the structure of 
the controlled unitary operators is shown in Fig. 8. The 
twofold controlled unitary operators in them are incom- 
pletely specified, and can be simplified by just turning off 
the operation of the first control qubit as shown in Fig. 8. 
Then all the twofold controlled gates in the circuit can 
be simplified. 



V. SCALING OF THE ALGORITHM 



composed of three gates: C\ C, and NOT gates) on each 
line of the circuit, except the last line, which has one gate. 
Therefore the total number of gates is N{n, 1) = 4n — 3. 
The number CNOT gates can be derived as follows: de- 
note A''cnot('t-, 1) as the number of CNOT gates. From 
Fig. 7, we have A^cnot('^, 1) - ^CNOT(n- - 1, 1) = 2(by 
adding one control qubit, one obtains two more con- 
trolled gates, a CNOT gate and a C-H gate). Then 
keep this decomposition procedure until n = 2, we ob- 
tain A'^cNOT(f^, 1) = 2n — 3. This is the upper bound of 
the number of CNOT gates for preparing a state in the 
one-electron system, since a target state in the incom- 
plete configuration space can be reduced to a state in the 
complete configuration space with fewer qubits. 

For the two-electron system, denote N{n, 2) as the 
number of gates needed to prepare the target state 
|4't(?t-, 2))from the jO)*^". Apply the decomposion pro- 
cedure in Fig. 6 until n = 2, and using the result from 
the one-electron system that N{n, 1) — An ~ 3, we ob- 
tain N{n, 2) 4n2 - lOn + 6. The number of CNOT 
gates can be derived as follows: let A^cnot('^, 2) repre- 
sent the number of CNOT gates needed to prepare the 
target state |'I't('^, 2)). Apply the procedure in Fig. 6 
until n = 2, and using the result from the one-electron 
system, we obtain A'^cnot('t-, 2) ~ 2n^ — 6n + 4. 

The scaling for preparing a general state of the m- 
electron system is derived as follows: denote iV(n, to) as 
the maximum number of gates, including both CNOT 
and single-qubit gates, needed to prepare the state of 
given values of n and to. We should stress that all two- 
fold controlled gates can be simplified to one-fold con- 
trolled gates (see Sec. IV). Looking at Fig. 6 and assum- 
ing that each block in it contains the maximum number 
of gates possible for the size of the block, we have: 

N{n, m) = 2N{n - 1, to - 1) + N{n - 1, to) + 2, (34) 

The first order derivative of N{n, to) with respect to n, 
where N and n are integers and the derivative is taken 
accordingly, is: 

dN{n,m) _ N{n,m) - N{n - l,m) 
dn n — (n — 1) 

= 2N{n-l,m-l) + 2. (35) 



In this section, we analyze the scaling of the algorithm 
for state preparation. The cost of a quantum algorithm is 
usually given by the number of CNOT gates used. Here 
we count all the twofold controlled gates as twoqubit 
gates since all twofold controlled gates can be simplified 
to onefold controlled gates, as shown in Sec. YE. 

In the one-electron system, for states in the complete 
configuration space, the scaling of the algorithm is de- 
rived as follows: let N{n, 1) denote the total number 
of gates, including both CNOT and single-qubit gates, 
needed to prepare the target state |^t("-, 1)) from the 
initial state jO)*^". From the quantum circuit shown in 
Fig. 7, one can see that there are four gates (a H gate is 



Using Eq. (34), the second order derivative of N{n,m) 



d'^N{n,m) d f dN{n,m) 



dn^ 



dn 



dn 



= — (27V(71-1,TO-1)) 

dn 

= 2^N{n-2,m-2) + 2^, (36) 

After (to — 1) steps of taking the derivative of N{n,m), 
we have: 

jV(n to) ^ _ + 1 1) + 2"-\ (37) 
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U{p-l,q) 



H - |0) 
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|0)®(/'-2) 



' (p - 1, g - 1) _ f/ ' (p - 1, 9 - 1)"' . (p - 1. "?) 



// [ - |o) 
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FIG. 8; Quantum circuit for two controlled unitary operators C-U{p,q) and C-U{p,q — 1) and their simplified versions. The 
two unitary operators in the top circuit are expanded in the second circuit, which is simplified in the bottom circuit. 



Using the scaling law given above for the one-electron 
case, i.e. N{ti — 7ti + 1, 1) = An — 4m + 1, we find that: 



d"W(n,m) 



X 4 : 



tm+l 



(38) 



We can now obtain the scaling behavior for preparing a 
general pure state of the m-electron, n-spin-orbital sys- 
tem by integrating Eq. (38) m times: 



Nin, m) cx ^—n^' 

m! 



(39) 



which scales polynomially with n, the number of qubits 
in the system. Since N{n, m) is the maximum number 
of gates needed to prepare a state with given values of n 
and TO, the scaling in Eq. (39) can be thought of as the 
upper bound of the number of gates needed to prepare 
the target state |^^'T(?^, tj)) from |0)®". The number of 
C-H gates is the majority of the circuit, e.g., as shown 
in Fig. 8; a C-H gate is composed of two single-qubit 
gates, and C, and a CNOT gate. The number of 
CNOT gates is comparable with the number single-qubit 
gates, both scaling as 0(2™n™/TO!). 

We now compare the scaling 0(2™n™/m!) with 0(2"), 
the number of gates needed to prepare an arbitrary 
n-qubit state. We look for the condition that makes 
2'"n™/m! smaller than 2". Using the Sterling approx- 
imation we have: 



(2nY 



2en 

TO 



(40) 



Then the following equations need to be satisfied: 

(^)'"<2", (41) 



TOl0g2 



2en 
m 



< n, 



TO log2 n < n. 



(42) 



(43) 



One can see that if to < n/log2?T., our algorithm is 
more efficient than the algorithms mentioned in Refs. [iGj 
and Then we compare the scaling of our algorithm 
with another algorithm introduced by Ortiz et al. [20| . 

0{N^n^). Assuming N ^ to ) ' looking for the 

condition that makes 2™n™/ml smaller than N'^n'^, we 
have: 



< 



to! 



< 



n 



(-)' 
Vto/ 



(44) 



(45) 



One can see that as long as n > 2m, our algorithm is 
more efficient than the algorithm in Ref . [20| . In quantum 
chemistry, usually n 3> m, i.e., many more orbitals are 
needed to describe a given number of electrons. The 
number of orbitals needed to describe a fixed number of 
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electrons depends on the accuracy of the calculation and 
the specific states that are investigated. For example, 
to study the highly excited states of a molecular system, 
more orbitals are needed than in study of the ground 
state. The number of orbitals needed and the number 
of electrons are therefore independent. So our algorithm 
gives an efficient way of preparing a general state for 
accurate simulation of a wide range of molecular systems. 



VI. APPLICATION TO SOME MOLECULAR 
SYSTEMS 

In this section, we apply the algorithm to two molec- 
ular systems: the hydrogen molecule and the water 
molecule. We prepare the multi-configurational self- 
consistent field (MCSCF) wave function of these 
two molecules. The MCSCF wave function is a linear 
combination of a number of CSFs, where each CSF is a 
symmetry-adapted linear combination of electron config- 
urations. 

The state function is a CSF, i.e., it is an eigen-function 
of the operators and 5^, where L and S are the or- 
bital angular momentum operator and the spin opera- 
tor, respectively. In quantum chemistry, one is usually 
most interested in the spin-states, i.e., the states that 
have certain spin multiplicity, of the molecular system. 
For molecules that have space symmetry, the electronic 
states of the molecule can be categorized into difii'erent ir- 
reducible representations of their point group. One only 
needs to perform calculations on states that belong to a 
certain irreducible representation. Considering these ef- 
fects, the number of configurations in the state function 
can be drastically below the size of the Hilbert space. 
The state function is spanned only in a subspace of di- 
mension much smaller than ^ ^ ^ • Therefore many con- 
figurations do not appear in the decomposition procedure 
shown in Eq. (24), where the state function is decom- 
posed into a |0) branch and a |1) branch. As a result, 
fewer gates are needed to prepare the state than the con- 
servative estimate given in Sec. V. Note that the algo- 
rithm itself is not simplified in any way when considering 
the symmetries of the molecule. 



A. Hydrogen molecule 



We use the MCSCF method with the cc-pVDZ basis 



set |31| to study the electronic structure of the hydrogen 
molecule. This gives 20 spin orbitals. Two electrons are 
distributed on these orbitals with the restriction of the 
spin multiplicity of the spin state. Considering the Daoh 
symmetry of the molecule, the ground state of H2 is ^S]+ 
is the Mulliken symbol that is used to label many- 
electron states; S, g, and the superscript "1" repre- 
sents angular momentum, parity, symmetry with respect 



to a vertical mirror plane perpendicular to the principal 
axis, and the spin multiplicity). 

For 20 qubits, the dimension of the whole Hilbert 
space is 2^" ~ 10^; for a fixed number of electrons, 
the two-electron, 20-spin-orbital system is of dimension 
20 \ 

2 1 = 190; considering the spin multiplicity and the 

space symmetry, the ground state MCSCF wave function 
is composed of at most 16 electron configurations. We 
can see that by considering the spin multiplicity and the 
space symmetry, the state space is drastically reduced. 

f 20 

To prepare the state function in the full space of ( ^ 

using our state preparation algorithm, we need roughly 
800 CNOT gates. To prepare the MCSCF wave func- 
tion that is composed of 16 electron configurations, by 
performing a numerical calculation, we found that only 
37 CNOT gates and 31 single-qubit gates are needed in 
total. This number, 37, is well below the conservative 
estimate of about 800 CNOT gates needed to prepare an 
arbitrary two-electron, 20-spin-orbital state. 



B. Water molecule 

For another example, the water molecule H2O, we use 
the MCSCF method with the cc-pVDZ basis set [3l|. For 
the ground state, considering the C2V symmetry of the 
water molecule, the Hartrec-Fock wave function of the 
H2O molecule is: 



(lai)2(2ai)2(162)2(3ai)2(16i)^ 



(46) 



The ground state of II2O is the ^Ai state. We apply 
a complete active space (CAS) type MCSCF method in 
order to reduce the cost of the calculation: the first two 
ai orbitals are frozen, the active space consists of the 
3ai-6ai orbitals, I61, I62 and 262 orbitals. So there are 
14 spin orbitals and 6 electrons in the active space. 

For 14 qubits, the dimension of the whole Hilbert space 
is 2^^ ~ 10^; the six-electron, 14-spin-orbital system is of 
14 



dimension 



6 



3003. To prepare the state func- 



tion in this space using our algorithm, we need roughly 
670000 CNOT gates, which is larger than 2" ~ 10"*, 
the number of CNOT gates needed to prepare a state in 
the full Hilbert space of 14 qubits. Considering the spin 
multiplicity and the space symmetry of the molecule, the 
ground state MCSCF wave function is composed of 152 
electron configurations. To prepare the state function us- 
ing our state preparation algorithm, performing the nu- 
merical calculation, we found that we need 1472 CNOT 
gates and 1146 single qubit gates. This represents a con- 
siderable and remarkable reduction of over 3 orders of 
magnitude in the number of CNOT gates (from ~ 670000 
to ~ 1470) needed to prepare an arbitrary six-electron, 
14-spin-orbital state. 

From these two examples, we can see that if the num- 
ber of electrons is comparable to the number of spin or- 
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bitals, the number of CNOT gates needed to prepare the 
corresponding state increases very fast. The algorithm 
is no longer efficient. However, in quantum chemistry, 
the number of orbitals needed for simulation purposes is 
usually much larger than the number of electrons. In the 
case of n ^ m, our algorithm provides an efficient way in 
preparing states of molecular systems, as demonstrated 
in the example of the hydrogen molecule above. 



VII. CONCLUSION 

In this paper, we present an efficient quantum algo- 
rithm for preparing a pure molecular-system-like state. 
The simulated system lies in the Fock space of a given 
number m of electrons and n spin orbitals on a quantum 
computer. The state wave function is a configuration 
state function and is entangled in general. In general 
this is to prepare a pure state in the combinatorial space 

^ ^ after the Jordan- Wigner transformation. 

In our algorithm, instead of starting from an initial 
state and transforming it to the target state, we start 
from the target state and transform it to the initial state. 
Then by inverting the quantum circuit sequence, we ob- 
tain the circuit for preparing the target state. A re- 
cursive procedure is employed for solving this problem. 
The twofold controlled gates that appear in each step 
of the recursive procedure can be simplified to one-fold 
controlled gates by turning off the operation of the first 
control qubit, since these gates are incompletely spec- 
ified. We show that at most 0(2'"+^n"'/m!) gates, in- 



cluding both CNOT and single-qubit gates, are needed to 
prepare a general state of the m-electron, n-spin-orbital 
system, which scales polynomially with n. The number 
of CNOT gates scales as 0(2™^™/?^!). This state prepa- 
ration algorithm works for systems of arbitrary number 
of electrons, but the number of CNOT gates needed will 
increase exponentially if the number of electrons m is 
proportional to the number of spin orbitals n. 

As examples, we have simulated our state preparation 
algorithm for the hydrogen and water molecules. In these 
two specific cases we analyzed, using the known sym- 
metries of the molecules, we found that the number of 
CNOT gates is reduced by up to 3 orders of magnitude. 
This provides a remarkable simplification to this type of 
problems. 

Note added. — In the final stages of preparing this pa- 
per, we became aware of a paper that outlined a quan- 
tum algorithm for the preparation of many-particle states 
on a lattice [13, in which the state is prepared in the 
first-quantized representation, while ours is in the second- 
quantized representation. We give an explicit and elabo- 
rate solution to the problem and in a very different way. 
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